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Abstract 

The gluon collision process that creates a heavy-quark-antiquark pair with small relative mo- 
mentum and large transverse momentum predicts at leading-order in the QCD coupling constant 
that the transverse polarization of the pair should increase with its transverse momentum. Mea- 
surements at the Fermilab Tevatron of the polarization of charmonium and bottomonium states 
with respect to a particular spin-quantization axis are inconsistent with this prediction. How- 
ever the predicted rate of approach to complete transverse polarization depends on the choice of 
spin-quantization axis. We introduce axes that maximize and minimize the transverse polarization 
from the leading-order gluon collision process. They are determined by the direction of the jet that 
provides most of the balancing transverse momentum. 

PACS numbers: 13.85.-t, 13.88.+e,13.90.+i 
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Important information about the fundamental interactions between elementary particles 
is encoded in the spins of particles produced in high-energy collisions. It is difficult to access 
that information, because the spin cannot be measured directly. The measurement of the 
spin state of a stable particle is usually not possible in a high-energy collider experiment. 
Information about the spin state of an unstable particle can be inferred from the angular 
distribution of its decay products. That angular distribution can be strongly correlated with 
the momenta of other particles in the final state, and if they are not measured, information 
about the spin is diluted. Hadron collisions have the additional complication that the 
fundamental interactions involve collisions of partons with longitudinal momenta that must 
be integrated over, which further dilutes the information about the spin. 

The accessible information about the spin of the unstable particle could in principle be 
obtained by measuring the complete angular distribution of its decay products for all possible 
values of the other relevant final-state variables. In practice, there is usually not enough data 
to make this feasible. The challenge is to find simple measurements that extract as much 
information as possible about the spin. In a two-body decay, the simplest measurement 
is the distribution of a polar angle 9 with respect to a specified axis. If that axis is used 
as the spin- quantization axis (SQA) for the unstable particle, the 9 distribution gives the 
probabilities for the projection of the spin along that axis. Additional information about 
the spin could be obtained by also measuring the polar angle distribution about a second 
SQA. These SQA's can be chosen as functions of other final-state variables to maximize the 
information about the spin. 

An interesting unsolved problem in high-energy physics is the spin dependence of the 
production rates for heavy quarkonia, whose constituents are a heavy (charm or bottom) 
quark and its antiquark. The nonrelativistic QCD (NRQCD) factorization framework ex- 
ploits the large mass of the heavy quark to express a cross section for inclusive quarkonium 
production as the sum of products of parton cross sections that can be calculated using per- 
turbative QCD and NRQCD matrix elements that can be measured in other experiments 
This approach leads to unambiguous predictions for the spin dependence of the cross 
sections. Measurements of the spin dependence of charmonium and bottomonium states 
produced at the Tevatron pp collider J2|j_ 3 Jj] are inconsistent with predictions based on 



leading-order parton cross sections 



8j . The resolution of the discrepancy may lie in 



large perturbative corrections to the parton cross sections 
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Most theoretical predictions and the experimental measurements of quarkonium polariza- 
tion have been carried out for a single SQA. In this paper, we point out that more detailed 
information about the spin can be extracted by measuring the polarization for a pair of 
SQA's for which the predicted polarizations are as different as possible. We construct a pair 
of optimal SQA's for spin-triplet S'-wave quarkonium states with large transverse momentum 
under the assumption that the dominant parton process is a gluon-gluon collision. 

To be specific, we consider the production of the spin-1 charmonium state J /ip in hadron 
collisions. There is a simple argument that a direct J/V', i-e. one that does not come from 
the decay of a heavier particle, should be increasingly transversely polarized as its trans- 
verse momentum Qt increases [11!]. Of the leading-order parton processes that produce a cc 
pair with large Qt, the gluon collision process gg — * {cc)g has the largest cross section. At 
asymptotically large Qt-, this cross section is dominated by gluon fragmentation, the produc- 
tion of an almost on-shell gluon by gg gg followed by the splitting of that virtual gluon 
into a coUinear cc pair. The cc pair is created in a color-octet state and is predominantly 
transversely polarized for almost any choice of SQA. Heavy-quark spin symmetry implies 
that the binding of this cc pair into a. J/ip is unlikely to change the spin states of the heavy 
quarks. Thus, the J/ ip should be increasingly transversely polarized as Qt increases. These 
qualitative arguments are supported by quantitative calculations using the NRQCD factor- 
ization formulas with leading-order parton cross sections . However, measurements 
by the CDF Collaboration seem to be completely incompatible with these predictions 

The longitudinal polarization 4- vector for a. J/ip with 4- momentum Q must satisfy 
Q ■ eL = and e| = —1. The most general 4- vector satisfying these conditions can be 
written in the form 



e^ = X7V-X2, X^ = {-g^'' + Q^Q''/Q')X,, (1) 

where X is a 4- vector. The physical interpretation of X is that in a J/-?/' rest frame, the 
direction of the 3- vector —X is the SQA of the J/ip. If X is timelike, the projection of the 
spin along the SQA coincides with the helicity in the rest frame of X. If the J ftp is produced 
in the collision of two hadrons with 4-momenta Pi and P2, the SQA is generally chosen to 
lie in the production plane defined by the momenta of the colliding hadrons in the J/ip rest 
frame. Thus, X in Eq. ([T]) has the form 

X" = aP^ + bP^, (2) 



where a and b are scalar functions. A simple example is the cm. helicity axis specified by 
-^cmh = -Pi + P2- The projection of the spin of the J/ip along this axis is its helicity in 
the center-of-momentum (cm.) frame of the colliding hadrons. In the J/ip rest frame, the 
SQA specified by Eq. ([2]) is antiparallel to the unit vector X = {aPi + bP2)/{\aPi + bP2\), 
so it is determined by the ratio a/b. If the J/ip is the only particle in the final state 
whose momentum is measured, then a/b can only depend on the 4- vectors Pi, P2, and Q. 
If additional information about the final state is measured, a/b can also depend on this 
information. 

If a cross section is dominated by a specific parton process, a natural prescription for 
a pair of optimal SQA's is that they maximize and minimize the transverse cross section 
from that parton subprocess at leading order in the QCD coupling constant a^. As a simple 
illustration of a pair of optimal SQA's, we consider dilepton production by the Drell-Yan 
mechanism qq in the parton model with intrinsic transverse momentum [l2, 



which can be regarded as a model for the effects of soft gluon radiation from the colliding 
partons. The cross section for a dilepton with longitudinal polarization vector given by 
Eqs. (HD and (El) is 



where Xi and X2 are the longitudinal momentum fractions of the colliding partons, is 
the electric charge of the quark, and {k"^) is the mean-square transverse momentum of the 
colliding partons. The cross section has been expanded to second order in the intrinsic 
transverse momentum vectors and then averaged over them. The longitudinal cross section 
in Eq. ([3]) depends on the ratio a/b. It is minimized by choosing a/b = —Xi/x2- It can 
be maximized by choosing a/b = +Xi/x2, so that the denominator vanishes. The values 
a/b = ~^X\/x2 that maximize and minimize the transverse cross section depend on the 
parton momentum fractions only through the ratio X\/x2- Under the assumption that the 
cross section is dominated by qq — > /i^/i^, we can derive an expression for X\/x2 in terms 
of variables that can be directly measured. At leading order in the intrinsic transverse 
momentum vectors, the energy-momentum conservation condition Q = XiPi + X2P2 implies 
Xi/x2 = Q ■ P2/Q ■ Pi- The 4- vectors X in Eq. ([2]) associated with the maximal and minimal 
SQA's can then be expressed as 





(3) 



The upper and lower signs correspond to the Collins-Soper axis 12|] and the perpendicular 



helicity axis introduced in Ref . , respectively. The perpendicular helicity axis is so named 
because the projection of the total spin of the dilepton along this SQA coincides with its 
helicity in the frame obtained from the hadron cm. frame by a longitudinal boost that makes 
the dilepton momentum perpendicular to the beam direction. Lam and Tung pointed out 
that the Collins-Soper axis maximizes the transverse polarization in the parton model with 



intrinsic transverse momentum 13|. That the perpendicular helicity axis minimizes the 
transverse polarization in this model does not seem to have been pointed out previously. 

We now consider the direct production of a. J/ip with large transverse momentum Qt ^ 
Aqcd- Since the mass of the charm quark is large compared to the scale Aqcd of 
nonperturbative effects in QCD, the inclusive cross section can be calculated using NRQCD 
factorization formulas [ll]. The leading-order parton processes that create a cc pair with 
small relative momentum are qq — » {cc)g, qg {cc)q, qg {cc)q, and gg — » {cc)g. The 
parton process with the largest cross section is gg {cc)g. This process is further enhanced 
by the growth of the gluon distribution at small values of the parton momentum fraction 
X as the momentum scale increases. We therefore focus on the parton process gg — * {cc)g. 
The color-octet ^5*1 state of the cc pair [cc8(^S'i)] dominates at asymptotic Qt, because it 
has a fragmentation contribution. At leading order in a^, the dependence of the longitudinal 
differential cross section on the scalar functions a and h defined by Eq. ([2]) is 

ddi a^xlu^ + fo^xft^ + (ax2 — hxiYs^ . . 

'^"rfSQ °^ (a«;i + bw2Y- abQ^s ' ^ 

where wi = Q ■ Pi, W2 = Q ■ P2, s is the cm. energy of the colliding hadrons, and s = X1X2S, 
t = 4m^ — 2xiWi, and u = 4m^ — 2x2^^2 are the parton Mandelstam variables. There is a 
delta function constraint on these variables: 

2{xiWi + X2W2) = X1X2S + Ami. (6) 

Minimizing and maximizing Eq. ([5]) with respect to a/b, we get 

Xi{s — i) 



33,max X2{s — u) ' 



(7) 



gg,mm X2[P (t — u) + U^{t — s)] 
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At asymptotic transverse momentum Qt ^ 2mc, the corresponding 4-vectors reduce to 

X^' + f9l 

^r.,min ^ ^ - ^ = ^CS- (10) 

Note that the minimal gg axis reduces in this hmit to the Colhns-Soper axis. 

Our optimahty criteria were based on the assumption that the parton process gg — *■ 
cc^i^Si) + g dominates. If that parton process imphes values for xi and X2 that are deter- 
mined by a measurable property of the final state, we can insert those values into Eqs. ([7]) 
and ([8]) to obtain optimal SQA's that are experimentally useful. The large transverse mo- 
mentum Qt of the cc pair is balanced by that of the recoiling gluon, which produces a jet of 
hadrons with nearly coUinear momenta. The polar angle of the recoiling gluon in the hadron 
cm. frame is approximately equal to the polar angle ^jet of the jet. The ratio Xi/x2 can be 
expressed as a function of ^^jct and the transverse and longitudinal momenta Qt and Ql oi 
the J /if) in the hadron cm. frame: 

xi _ {Qo + Ql) singjet + Qt(1 + cos^jet) 
X2 {Qo - Ql) sin6'jet + Qt{1 - cos6'jet) ' 

where Qo = (Qt + Q\ + ^j/^,)^^^- Given the value of Xi/x2, the delta function constraint 
in Eq. ([6]) can be solved for xi and X2- The two-fold ambiguity associated with the root of a 
quadratic polynomial is resolved by choosing the smaller values of Xi and X2- Inserting them 
into Eqs. ([7]) and ([8]), we obtain expressions for a/h that depend on quantities that can be 
directly measured. We will refer to the corresponding SQA's as the maximal and minimal 
gg axes, respectively. The prescriptions for these optimal SQA's can be extended beyond 
leading order in as by choosing ^jet in Eq. (|TT1) to be the angle in the hadron cm. frame 
of the jet with the largest transverse energy. Since the direction of the jet is insensitive to 
soft gluon radiation and to the splitting of a parton into collinear partons, QCD radiative 
corrections to the polar angle distributions defined by our optimal SQA's can be calculated 
systematically using perturbative QCD. 

To illustrate the use of pairs of SQA's, we consider direct J /ip production at the Tevatron, 
which is a pp collider with cm. energy 1.96 TeV. We use NRQCD factorization formulas with 
leading-order parton cross sections to calculate the transverse and longitudinal cross sections 
for J /if). For the NRQCD matrix elements, we use central CTEQ5L values from Ref. [3] with 
X = 1/2: {Oi'^{^Si)) = 1.4 GeV^ {Ol'^{^Si)) = 0.0039 GeV^ {Oi'^{^So)) = 0.033 GeV^ 




FIG. 1: Polarization variable a for various spin-quantization axes as functions of Qt/Mji^ for 
direct J/ifj at the Tevatron (left panel) and direct T(IS') at the LHC (right panel). The SQA's 
are the maximal and minimal gg axes (upper and lower solid lines), the perpendicular helicity and 
Collins-Soper axes (upper and lower dashed lines at large Qt)-, and the cm. helicity axis (dotted 
line). 

and {Og^^{^Po)) /ml = 0.0097 GeV^. In the parton cross sections, we set rric = Mj/^/2 = 
1.5 GeV. We use the CTEQ6L parton distributions [15] with 3 flavors of quarks and the 
next-to-leading-order formula for as{fi) with 4 flavors of quarks and Aqcd=326 MeV. The 
factorization and renormalization scales are set to fi = {Q^ + Q^V^^- We impose a rapidity 
cut \y\ < 1 on the J/ip momentum. 

A convenient polarization variable for J/ip is a = {(Tt — 2aL)/{aT + 2crj^), whose range 
is — 1 < a < +1. We consider flve SQA's: the maximal and minimal gg axes deflned by 
Eqs. ([7]) and ([H]) together with Eqs. (E]) and ( fTTl) . the perpendicular helicity and Collins- 
Soper axes deflned by X^h and Xcs in Eq. (jlj), and the cm. helicity axis, which has been 
used in most previous work on this problem. The leading-order predictions for a for the flve 
SQA's are shown as functions of Qt/Mji^ in the left panel of Fig. [H For the cm. helicity 
and perpendicular helicity axes, a increases with Qt and asymptotically approaches 1. The 
two axes are essentially identical at large transverse momentum. For the maximal gg axis, it 
approaches 1 much more rapidly. For the Collins-Soper axis, a decreases with Qt, reaching 
—0.33 at Qt = 10 Mjj^. For the minimal gg axis, a approaches the same asymptotic value 
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much more rapidly. 

We also consider direct T(15') production at the CERN LHC, which is a pp collider with 
cm. energy 14 TeV. For the NRQCD matrix elements, we use central CTEQ5L values from 
Ref. 0: {Oj(^Si)) = 10.9 GeV^ {Oj{^Si)) = 0.025 GeV^ {Oj{^So)) = 0.068 GeV^ and 
{Oj{^Po))/ml = 0.014 GeV^. We have used the averages of the values for the (Oj(^S'o)) = 



and {Oj(^Po)) = cases of Ref. 



16|. In the parton cross sections, we set rrib = Mr/2 



4.7 GeV. We imposed a rapidity cut \y\ < 3 on the T. The leading-order predictions for a 
for the five SQA's are shown as functions of Qt/M-^ in the right panel of Fig. [H 

Our optimal SQA's are not useful for fixed-target experiments, because the recoiling 
jet is usually not observed. However, the angle ^jct of the recoiling jet can be measured 
relatively easily in high-energy hadron colliders. To provide some idea of how much the 
data sample will be decreased by the requirement that ^jet be measured, we impose a cut 
on the pseudorapidity r^jet = lntan(6'jet/2) of the jet. For J/ip at the Tevatron, the fraction 
of events satisfying the J/ip rapidity cut \y\ < 1 that also survive a jet pseudorapidity cut 
|?7jet| < 1 is greater than 0.21 for Qt > Mj/^. For T(15') at the LHC, the fraction of events 
satisfying the T rapidity cut that also survive a jet pseudorapidity cut |?7jet| < 3 is greater 
than 0.72 for Qx > M^. These fractions are large enough that measuring 6'jet should not 
dramatically decrease the size of the data sample. 

The CDF and DO Collaborations have measured the polarization as a function of Qt for 
charmonium and bottomonium mesons produced at the Tevatron For J/ip, the 

variable a for the cm. helicity axis was measured for Qt/Mj/^ as high as 9.7. For T(IS'), 
a for the cm. helicity axis was measured for Qt/My as high as 2.1. At the LHC, it should 
be possible to measure a for these and other heavy quarkonium mesons out to much larger 
values of Qt- 

The measurement of a for two different SQA's will provide more information about 
the spin if there is a large difference in the prediction of the transverse polarization with 
respect to the two axes. Our leading-order calculations predict a large difference in a for the 
perpendicular helicity and Collins-Soper axes for Qt ^ M. One advantage of these two axes 
is that Xcs and in Eq. (j4]) are invariant under independent longitudinal boosts of the 
two colliding hadrons. This implies that a for these axes is insensitive to coUinear radiation 
from the colliding partons, so radiative corrections may be smaller than for the cm. helicity 
axis. Our leading-order calculations predict an even larger difference in a for the maximal 
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and minimal gg axes. These optimal axes require the measurement of the direction of a 
recoiling jet, but we have shown that this should not dramatically decrease the size of the 
data sample. 

Similar methods could be used to derive optimal SQA's for the production of heavy el- 
ementary particles in the standard model, such as the weak bosons and and top 
quark. There should be a large difference in the polarization between the perpendicular 
helicity and CoUins-Soper axes, but an even larger difference between the appropriate max- 
imal and minimal SQA's. If predictions for the polarization of these particles with respect 
to optimal SQA's can be verified, then optimal SQA's will also provide a new window into 
the spins of new particles created at the LHC. 
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